Occupation Statistics of a BEC for a Driven Landau-Zener Crossing 
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We consider an atomic Bose-Einstein condensate (BEC) loaded in a biased double-well trap with 
tunneling rate K and interatomic interaction U. The BEC is prepared such that all N atoms are 
in the left well. We drive the system by sweeping the potential difference £ between the two wells. 
Depending on the interaction u = NU/K and the sweep rate £, we distinguish three dynamical 
regimes: adiabatic, diabatic, and sudden and consider the occupation statistics of the final state. 
The analysis goes beyond mean-field theory and is complemented by a semiclassical picture. 



The theoretical and experimental study of driven 
atomic Bose-Einstein Condensates (BEC) in a few site 
system using optical lattice technology has intensified in 
recent years [J 0, [3, 0]- Beyond the fundamental in- 
terest of these studies, they also aim to create a new 
generation of nanoscale devices such as atom transistors 
[5(. Consequently, a wealth of research has been done 
on the prototype two-site (dimer) system, either within 
the framework of a nonlinear mean- field approach [l| , op- 
tionally using higher order cumulants la , or adopting a 
conventional many-body perspective @, 3, 0, [l(| ■ Such 
investigations have revealed many interesting phenom- 
ena associated with eigenvalue spectra, the structure of 
the eigenstates, wave-packet dynamics, e.g. of the Bloch- 
Josephson type, and leaking dynamics due to dissipative 
edges. 

Driven dime rs p rove to be even more challenging 



ill lia lid Uilia • 

The scenario investigated in these 



studies involves many-body Landau-Zcner (LZ) transi- 
tions induced by sweeping the potential difference £ be- 
tween the two wells. Specifically, one assumes that ini- 
tially £ is very negative and that all n — N atoms are in 
the first well. Then £ is increased at some constant rate 
£ to & very positive value. The objective is to calculate 
how many atoms (n) remain in the first well. The major- 
ity of published works are based on the Gross-Pitaevskii 
equation (or its discrete analogue), which is a mean- field 
approach 11, [e|. There are only a few studies that 
have made further progress within the framework of a full 
quantum mechanical treatment of the system 1^, 14 , 15| . 
However, they all focus on calculating the average occu- 
pation (n) , neglecting to form a theory for the occupation 
statistics P(n) and, in particular, for the variance Var(n). 

Outline. — In this Letter, we consider a driven dimer 
with intersite hopping amplitude K and interatomic in- 
teraction U. The parameter u = NU/K can be either 
positive (repulsive) or negative (attractive) and its mag- 
nitude determines various dynamical regimes Rabi 
(|u| < 1), Josephson (1 < |u| < iV 2 ), or Fock > N 2 ). 
Depending on the sweep rate £ we distinguish between 
adiabatic, diabatic, and sudden dynamical scenarios and 
study the asymptotic occupation statistics as a function 



of u. Our analysis goes beyond mean- field theory, using a 
semiclassical picture and involving detailed simulations. 

Modeling. — The simplest model that describes inter- 
acting bosons on a lattice is the Bose-Hubbard Hamilto- 
nian (BHH), which in case of the dimer (two sites) reads: 



H= 



2 r 

£ 

i=i 



£ihi + —n,i(hi - 1) 



K 



with H=l, where bi and b\ are bosonic annihilation and 
creation operators and hi — b\bi counts the number 
of particles at site i = 1,2. The validity of this two- 



mode ap pro ximation 171 ] for a double well is discussed 



in Refs. [18|, ll9|, as well as in Ref. [20| for biased sys- 



tems, and it was found to yield good agreement with 
the experiment (l6| . The total number of particles N = 
hi + hi is a constant of motion, allowing us to consider a 
Hilbert-space of dimension M = N+l, which is spanned 
by the Fock basis states \n) = \ni=n,n2=N—n). Below 
we assume an even N ^> 1 and define j = N/2, hence 
TV = 2j+l. The BHH for a given N is formally equiv- 
alent to the Hamiltonian of a spin j particle. Defining 
J z = (hi — 7i2)/2 and J+ = b\i)2 it can be rewritten as: 



H = U J z + £J Z — KJ X ; £=£i-£ 2 



(2) 



where £ is the bias. In the absence of interaction, this 
Hamiltonian can be reinterpreted as describing a spin in a 
magnetic field with precession frequency Q = (-K, 0, £). 

The many-body Landau-Zener scenario assumes that 
initially all particles are located in the first site 
|\I/(f=0)) = \N). The bias £ is then varied with some 
constant sweep rate £. At the end of the sweep the 
occupation statistics P{n) = |(n|\&(t))| 2 becomes time- 
independent. Depending on the outcome (Fig. 1) we dis- 
tinguish between: (i) an adiabatic process P(n) ~ 5 n ,o', 
(ii) a sudden process P(n) « <5n,Ar; and (iii) a diabatic 
process P(n) as 8 n ,n cl where n c ^ 0, N. 

Phase Space. — In order to analyze the dynamics for 
finite U, it is convenient to rewrite the BHH using canon- 
ical variables. Formally, our system corresponds to two 
coupled oscillators and thus we can define action-angle 
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FIG. 1: (Color online) Phase space and corresponding energy 
levels (conceptual plot). Note that a ground state preparation 
of a U < system with increasing bias £(t) is equivalent to a 
preparation in the most excited state of a U > system with 
decreasing bias. The latter convention has been adopted in 
the right panel and in our simulations. 



variables bi = yWi exp(zy)j). Note that the translation- 
like operator exp(itp) is in fact non-unitary because it 
annihilates the ground state, but this is irrelevant for 
N > 1 [HHH. With these coordinates the BHH takes 



a form that resembles the Josephson Hamiltonian: 



n 



NK 



u(cos 9) 2 + e cos 9 — sin 9 cos <p 



(3) 



where 9 is an alternative way to express the oc- 
cupation difference J z = (ni— n 2 )/2 = (iV/2)cos(0) and 
ip = ipi~ip2- The scaled bias is e = £/K. Note that (p 
and 8 do not commute. The classical phase space is de- 
scribed either using the canonical coordinates (ip, n), with 
n = J z +(N/2) G [0,7V], or equivalently using the spheri- 
cal coordinates (ip, 9) . In the former case the total area 
of phase space is 27riV with a Planck cell 2irH and h=l, 
while in the latter case the phase space has total area 
4-7T with a Planck cell Att/N . Within the semiclassical 
approximation, a quantum state is described as a dis- 
tribution in phase space and the eigenstates are associ- 
ated with stripes that are stretched along contour lines 
TC(<p, 9) — E. The energy levels E n can be determined 
via WKB quantization of the enclosed phase space area. 

Separatrix. — The phase space topology becomes 
non-trivial, i.e. has more than one component as in 
Fig. la, if |u| > 1 and |e| < e c where e c = (u 2/3 - 1) 3/2 . 
In this regime, a separatrix divides the phase space into 
three regions: two islands (green) that contain the up- 
per energy levels and a sea (blue) that contains the lower 
energy levels. In this description and in the numerics 
below we assume U > and adopt the following enumer- 
ation convention: we define Eq as the most excited level, 
while En corresponds to the lowest one. Note that the 
replacement U i— > —U would merely invert the order of 
the levels, and E would become the ground state. 

Fig. lb is an illustration of the energy levels as a func- 
tion of the bias. Each level is associated via WKB quan- 
tization with one of the contour lines in Fig. la. For the 
sake of our later analysis we define E nc as the level which 
is closest to the separatrix for |e| = e c . At this critical 
value of the bias one of the islands has a vanishing phase 
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FIG. 2: (Color online) Dynamical evolution in a situation 
in which the mean-field (Gaussian) approximation does not 
hold: iV=100 particles are prepared in one site of a sym- 
metric (£=0) double well and then evolved with a u=2 
Hamiltonian. The initial wave-packet is stretched along the 
separatrix. (a) (semi) classical distribution in phase space. 

(b) Wigner function of the corresponding quantum state. 

(c) Quantum P(n) compared with the (semi)classical and 
with the mean-field (binomial/Gaussian) predictions. The 
average occupation (n)=47 is associated here with a huge 
super- binomial variance Var(n)=1530 instead of the binomial 
value Var(n)^25. In all the simulations we use a 4th order 
Runge-Kutta and made sure that during the integration the 
probability leakage is <C 10 -6 . 



space area, while the area of the other is A c w A~Ke c /u 
Using WKB quantization we get 



4tt/N 



(1 - u- 2 ' 3 f/ 2 N, 



(4) 



where the approximation for A c has been derived using 
methods as in Ref. [Il[ , and has been tested numerically. 

Simulations. — As a preliminary step we did simu- 
lations of the undriven wave-packet dynamics with £ =0 
and u=2. For this value of u the separatrix crosses the 
north pole of the phase space and therefore the wave- 
packet stretches along it, as illustrated in Fig. 2a (clas- 
sical) and in Fig. 2b (quantum mechanical). We note 
that this type of dynamics cannot be properly addressed 
by the mean-field approximation. The mean-field equa- 
tion merely describes the Hamiltonian evolution of a sin- 
gle point in phase space and therefore assumes that the 
wave-packet looks like a minimal Gaussian at any mo- 
ment. Whenever the motion takes place near the sep- 
aratrix, the mean-field description becomes inapplicable 
and consequently the distribution P(n) is likely not to 
be binomial (Fig. 2c). In what follows we use the terms 
sub/super-binomial in order to refer to a P(n) with a 
smaller/larger spreading than the mean-field results. For 
the dynamics described in Fig. 2, the stretching along the 
separatrix leads to a super-binomial result. 

Next, we address the effect of separatrix motion on 
P(n) in the bias-sweep scenario. Note that this separa- 
trix motion cannot be avoided: For e < —e c the wave- 
packet is localized in the upper level. When e — —e c the 
separatrix emerges. As long as — e c < e < the wave- 
packet remains trapped in the top of the big island which 
gradually shrinks. When e becomes larger than zero, the 
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FIG. 3: (Color online) (a) Parametric evolution of the adia- 
batic energy levels versus the decreasing £i(t) for N=10 par- 
ticles and tt=2.5. We use units of time such that K=l. The 
dashed line corresponds to the level n c —4. (b) Average occu- 
pation (n) versus the sweep rate £. (c) Participation number 
(PN) and Var(n) versus £. The vertical lines indicate the 
various adiabatic and diabatic (dashed) thresholds. 



wave-packet can partially tunnel out from the shrink- 
ing island to the levels of the expanding island. When 
e = +s c , the shrinking island disappears and the remain- 
ing part of the wave-packet is squeezed out along the 
n = n c contour, resembling the dynamics of Fig. 2. One 
observes that the stretching along the separatrix during 
the nonlinear LZ transition is accompanied by narrow- 
ing in the transverse direction. This leads to a sub- 
binomial rather than super-binomial result for the dis- 
tribution P(n) at the end of the sweep. 

In Fig. 3 we plot the average occupation (n) and the 
participation number PN= [J2 n P(n) 2 ] -1 of the distribu- 
tion at the end of the sweep as a function of e. These, un- 
like (n), provides significantly more information regard- 
ing the nature of the crossing process. For very slow 
rates, the wave-packet follows a strict adiabatic process 
ending in n=0, i.e. all particles move to the other site. 
For a moderate sweep rate the wave-packet ends in a su- 
perposition of n=0 and n=l states, indicated by PN=2. 
We also resolve the possibility of ending entirely at n=l 
or at n=2 or at n=3. In the case shown in Fig. 3, we have 
n c rj 4. For larger sweep rates, we observe a qualitatively 
different behavior that can be described as a crossover 
from an adiabatic/diabatic behavior to a sudden behav- 
ior at the peak value PN=4. In order to appreciate the 
deviation of the numerical results from the mean-field 
theory prediction, we plot Var(n) versus (n) in Fig. 4 
and compare with the binomial expectation. We further 
analyze the observed results in the last section. 

Thresholds. — The various thresholds that are in- 
volved in the adiabatic-diabatic-sudden crossovers are in- 
dicated in Fig. 5. They all follow from the breakdown of 
a "slowness condition" that can be written as 
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FIG. 4: (Color online) Y = Var(n)/iV versus X = (n)/N. 
Symbols correspond to numerical data, while lines indicate 
the sub-binomial scaling relation Eq.© with u—2.5 (dotted) 
and w=4.05 (solid). No fitting is involved. As a reference we 
plot the standard binomial scaling (dashed line) which would 
strictly apply in the absence of interaction (u=0), and the 
results of semiclassical (SC) simulations which are obtained 
by running ensemble of trajectories in phase space. There is 
a clear crossover from a binomial to the sub-binomial scaling 
and the agreement becomes better for large N. Note that 
the super-binomial data point (X, Y) = (0.47, 15.3) that cor- 
responds to the distribution in Fig. 2 is out-of-scalc. 



where u> OBC is a characteristic frequency of the unper- 
turbed dynamics and k is the coupling parameter that 
determines the rate of the driven transitions. 

In the strict quantum adiabatic framework, w OBC is sim- 
ply the level spacing and k is determined by the slopes 
of the intersecting levels. In order to determine the adi- 
abatic thresholds in Fig. 3 we observe that for the in- 
tersection of the 0th level with the (N—n) level the dif- 
ference in slope is k = (N—n), because asymptotically 
d(E n —E m )/d£ ~ (n— to). In the absence of interaction 
(it <C 1), the level spacing is lu osc — K and only nearby 
levels are coupled, leading to the standard Landau-Zencr 
adiabaticity condition £ <C K 2 . With strong interaction 
there is an Nth order coupling between the n=0 level 
and the n=N level, which allows tunneling from the 
top of one island to the top of the other island (as il- 
lustrated in Fig. 1). An estimate for this coupling is 
K e{! = [NK\/^ N - 1 {N-1)\\{K/U) N - 1 0. Similar con- 
siderations can be applied to the bottom sea level leading 
to the distinction between mega, gradual, and sequential 
crossings (see Ref. 22]). 

For large N it might be practically impossible to sat- 
isfy the strict adiabatic condition which is associated with 
the possibility to tunnel from the top of one island to the 
top of the other. Then the relevant mechanism for tran- 
sition, i.e. the emission to the level n c as described in the 
previous section, becomes semiclassical. The frequency 
that governs this process is the oscillation frequency at 
the bottom of the sea a-i osc ~ \NU K^/ 2 . It determines 
the level spacing of the lower energy levels and also de- 
scribes the level spacing in the vicinity of the separatrix, 
apart from some logarithmic corrections [23]. It follows 
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FIG. 5: (Color online) Diagram of the (U, £) regimes for a 
ground state preparation. In the Rabi regime \U\ < K/N 
we have a crossover from adiabatic to sudden behavior. For 
K/N < U < NK (or U > NK) we have a broad crossover 
from adiabatic gradual (or sequential [z|) behavior to sud- 
den behavior. For U < —K/N we have two crossovers: the 
first from quantum adiabatic to diabatic behavior and the sec- 
ond from diabatic to sudden behavior. The diabatic behavior 
can be regarded as a classical nonlinear adiabatic behavior. 
For U < —NK the distinction between the diabatic and the 
sudden regime is blurred because the final state is the same. 

that the diabatic-sudden crossover involves the threshold 
condition £ -C \NUK\ as indicated in Fig. 5 and in Fig. 3 
for the specific parameters of the simulations. 

Scaling. — Due to the squeezing along the separa- 
trix, the spreading of the wavepacket for an idealized 
diabatic process becomes negligible in the transverse di- 
rection. The diabatic-sudden crossover is related to the 
non-adiabatic transitions between the remaining (N—n c ) 
sea levels, where nonlinear effects are negligible. It fol- 
lows that the spreading can be approximately modeled by 
the toy Hamiltonian TL — h{t) ■ J , where J is a spin entity 
with j off = (N—n c )/2, and h(t) is a field with constant 
magnitude \h(t)\ — u oso corresponding to the mean-level 
spacing. The sweep is like a rotation of h(t) in the plane 
with some angular rate oj. For such a (linear) model 
the mean-field approximation is exact and therefore we 
suggest (due to the truncation of Hilbert space) a sub- 
binomial rather than binomial scaling relation between 
the mean and the variance of the occupation statistics: 

Y = (1-X)|_i ; with C=^, (6) 

where X = (n)/N and Y = V&r(N)/N . Our numerical 
data is reported in Fig. 4 together with the binomial 
(c=0) and sub-binomial scaling relation Eq.((6]). The nu- 
merics confirm the expected u-dependent crossover from 
binomial to sub-binomial statistics, where the latter, with 
no fitting parameters, sets a lower bound for the variance. 

Summary. — In view of the strong research interest 
in counting statistics of electrons in mesoscopic devices, 
it is surprising that the issue of occupation statistics of 
BECs has been explored only for equilibrium phase tran- 
sitions. We were motivated to address this subject in the 
framework of dynamical processes by state-of-the-art ex- 
periments aimed at counting individual particles @, S Q • 



We have shown that in the case of a many-body Landau- 
Zener transition the mean-field binomial expectation is 
not realized, however a sub-binomial scaling relation still 
works quite well. The study of the occupation statistics 
and, in particular, the participation number of the final 
distribution, allowed us to resolve all the details of the 
adiabatic-diabatic-sudden crossovers and to verify theo- 
retical estimates for the threshold of each crossover. 

Acknowledgments. — This research was funded by a 
grant from the US-Israel Binational Science Foundation 
(BSF), the DFG Forschergruppe 760, and a DIP grant. 



[1] M. Albiez, et al., Phys. Rev. Lett. 95, 010402 (2005). 
[2] C.-S. Chuu, F. Schreck, TP. Meyrath, J.L. Hanssen, 
G.N. Price, and M.G. Raizen, Phys. Rev. Lett. 95, 260403 

(2005) ; A.M. Dudarev, M.G. Raizen, and Q. Niu, Phys. 
Rev. Lett. 98, 063001 (2007). 

[3] S. Foelling, S. Trotzky, P. Cheinet, M. Feld, R. Saers, A. 

Widera, T. Mueller, I. Bloch, Nature 448, 1029 (2007) 
[4] P. Cheinet, S. Trotzky, M. Feld, U. Schnorrberger, M. 

Moreno-Cardoner, S. Foelling, I. Bloch. larXiv:0804.3372l 
[5] J. A. Stickney, D. Z. Anderson, A. A. Zozulya, Phys. Rev. 

A 75, 013608 (2007). 
[6] J. R. Anglin and A. Vardi, Phys. Rev. A 64, 013605 

(2001); A. Vardi, V. A. Yurovsky, and J.R. Anglin, Phys. 

Rev. A 64, 063611 (2001); A. Vardi and J. R. Anglin, 

Phys. Rev. Lett. 86, 568 (2001). 
[7] G. Kalosakas, A. R. Bishop, and V. M. Kenkre, Phys. 

Rev. A 68, 023602 (2003); G Kalosakas, A R Bishop and 

V M Kenkre, J. Phys. B: At. Mol. Opt. Phys. 36, 3233 

(2003); G. Kalosakas and A. R. Bishop, Phys. Rev. A 65, 

043616 (2002). 

[8] F. Trimborn, D. Witthaut, H. J. Korsch, 

larXiv:0802.lT42l 
[9] M. Hiller, T. Kottos, and A. Ossipov, Phys. Rev. A 73, 

063625 (2006). 

[10] E. M. Graefe, H. J. Korsch, A. E. Niederle, Phys. Rev. 

Lett. 101, 150408 (2008). 
[11] J. Liu, L.-B. Fu, B.-Y. Ou, S.-G. Chen, and Q. Niu, Phys. 

Rev. A 66, 023404 (2002). 
[12] B. Wu and Q. Niu, Phys. Rev. A 61, 023402 (2000). 
[13] D. Witthaut, E. M. Graefe, and H. J. Korsch, Phys. Rev. 

A 73, 063609 (2006) . 

[14] P. Solinas, P. Ribeiro, R. Mosseri, larXiv:0807.0703l 
[15] A. Atland, V. Gurarie, Phys. Rev. Lett. 100, 063602 

(2008). 

[16] R. Gati and M. K. Oberthaler, J. Phys. B 40 (2007) R61. 
[17] G. J. Milburn, J. Corney, E. M. Wright, and D. F. Walls 

Phys. Rev. A 55, 4318 (1997). 
[18] R. W. Spekkens and J. E. Sipe, Phys. Rev. A 59, 3868 

(1999). 

[19] D. Ananikian and T. Bergeman, Phys. Rev. A 73, 013604 

(2006) . 

[20] D. R. Dounas-Frazer, A.M. Hermundstad, and L. D. 

Carr, Phys. Rev. Lett. 99, 200402 (2007) 
[21] A. J. Leggett, Rev. Mod. Phys. 73, 307 (2001). 
[22] M. Hiller, T. Kottos and D. Cohen, Europhys. Lett. 82, 

40006 (2008); M. Hiller, T. Kottos and D. Cohen, Phys. 

Rev. A 78, 013602 (2008). 
[23] E. Boukobza, M. Chuchem, D. Cohen and A. Vardi, 

Phys. Rev. Lett. 102, 180403 (2009). 



